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ABSTRACT 


To improve the stability and containment time of a plasma in a 
magnetic field, the 7 l?sma stabilization for controlled thermonuclear 
fusion must be developed. The methods ol principles of magnetic 
shear, magnetic well, dynamic, and feedback stabilization are of primary 
importance and will be given primary consideration, fhe objectives of 
this research are to: 

« review various stabilizing methods, 

« compare the effectiveness, advantages and disadvantages of 
some stabilizing methods, and 

• deduce from the results obtaine l a new idea or original 

method for understanding plasma instability and improving 


plasma stabilization. 




1. INTRODUCTION 


Plasma stabilization during developing controlled fusion is one of 
the most difficult achievements in thermonuclear research. To improve 
the stability and containment time of a plasma contained in a magnetic 
field, various stabilizing methods have been and are being developed. 

Among them, of primary importance, are magnetic field configuration 
stabilization, dynamic stabilization and feedback stabilization. For 
magnetic field configuration stabilization, the plasma may be made 
self-stabilizing in some complex magnetic fields. For dynamic stabili¬ 
zation, the high-frequency (or radiofrequency) field interacts with 
the entire plasma to establish a new dynamically stable equilibrium. 

For feedback stabilization, the plasma equilibrium remains the same 
and the time-varying suppression signal interacts only with the unstable 
part of the plasma. The reality of developing a controlled fusion reactor 
depends mainly on the development of an appropriate magnetic field con¬ 
figuration and the availability of remote control and detection methods 
of the dynamic stabilization or feedback stabilization for a Thermonuclear 
plasma. 

There are many methods proposed for the magnetic field configuration 
stabilization, the dynamic stabilization and the feedback stabilization 
of a thermonuclear plasma. Only some of the main methods of plasmc. stabil 
zation will be discussed, analyzed and compared in this research work. 
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2. MAGNETIC FIELD CONFIGURATION STABILIZATION 

Plasma stabilization by magnetic field configuration consists 
mainly of the following methods: (a) internal axial magnetic field; 

(b) rotational transform devices; (c) screw pinch; (d) magnetic shear; 

(e) magnetic well-minimum B (open-end system), magnetic well-minimum 
average B (clcsed-end system); and (f) multipole configuration. 

2.1 Internal Axial Magnetic Field 

To improve the stability of a linear or toroidal pinched plasma, 

1 2 

an internal axial magnetic field was early introduced. ’ Experiments 
show that the introduction of an axial magnetic field in plasma can delay 
instability of the pinched discharge. Certain ratios of the axial and 
azimuthal magnetic fields, b z /®q» and some magnitude of the pinch ratio, 
r Q /r, (less than 5) must also be met, where b ^ and Bg are the axial and 
azimuthal magnetic fields and r Q and r are the radii of the discharge 
tube and the pinched plasma, respectively. 

Interdiffusion (or intermixing) of axial and azimuthal magnetic 
fields must ba taken into consideration when B z ard B,^ diffuse into each 
other in the presence of a finite plasma resistivity, T)• The interdiffusion of 
the magnetic fields has two important effects:, (a) The redistribution of mag¬ 
netic field is in a transient state. Aj a result, the pinched plasma tends to 
become unstable, due to continuous decrease in B z< (b) The magnetic field 
diftui ion energy loss due to resistivity (1)) in the current sheath tends to heat 
the plasma, i.e., the ohmic heating. 
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Consider a cylindrical plasma sheath of unit length and radius r, 
contained by the axial and azimuthal magnetic fields in a long discharge 
tube. To evaluate the field time constant, t, for diffusion of the axial field 
flux out of the cylindrical sheath, the energy, W, of magnetic field B 
(in MKS units) is given by 


where ^ 
to time, 




is the permeability in vacuum. By differentiating B with respect 
t, the rate of the energy loss is expressed as 


oW _ _ nr 2 dB 2 
dt 2n o dt 


( 2 . 2 ) 


This relation is equal to resistive power, P, in the current sheath to 
heat the plasma 




(2.3) 


where I and R are the electrical current and resistance of the plasma 
and 6 is the sheath thickness. Assume that the rate of magnetic energy 
loss by diffusion is entirely converted into the resistive power of ohmic 
heating in the plasma, from Eqs. (2.2) and (2.3), it follows that 


2 _ _ r6 dB 2 
41)n o dt 


(2.4) 


A plasma sheath can be considered as a single-turn solenoid, according 
to the Biot-Savart law that current and magnetic field are related by 


H = yl 


(2.5) 
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Introduction of Eq. (2.5) into Eq. (2.4) gives 
. 2 p, r6 


_1_ dlT 

„2 at 


411 


Integrating this and using the initial condition B = B q at t = 0, the 

2 

time constant T and the field time constant T related to B and B are. 
o 

respectively 

-t/T 


B = B 


-t/T 


( 2 . 6 ) 


where 


p. r6 

T o * "4 tT ’ 


V 6 

211 


' that t = 2t 


It is obvious that the diffusion of the confining magnetic field (or the 
drift of charged particles) from the plasma can be slowed down if Tj is 
small and r is large. Tnis result can be generalized for various devices:’ 
(a) it will take a longer time for to diffuse in a large discharge tube 
than in a small one, and (b) the containment of a stable discharged plasma 
requires rapid initial heating to attain a high temperature for small Tj. 

2.2 Rotational Transform Devices 

It is well known that the rotational transform was first introduced 
in the stellarator machine.'* The stellarator can produce a rotational 
transform because of its "figure 8" configuration. In the presence of an 
axial magnetic field with a rotational transform, the electrostatic fields 
and charge separation due to the oprositeiy directed drifts of positive 
ions and negative electrons in the rtellarator system can he considerably 
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suppressed. Magnetic transform and particle transform can also serve 

4 

the purpose of the rotational transform. 

To compensate for the ion and electron drifts in the curved region 
of c race track type of toroidal device, a corrugated magnetic field can 
also produce the effect of a rotational transform. The corrugated mag¬ 
netic fields can be realized by a suitable spacing of the winding of 
a number of solenoidal coils which produce the axial magnetic field when 
energized. 

2.3 Screw Finch 

Screw pinch is a dynamic pinch wherein an external axial magnetic 
field, B z> is trapped outside the plasma. This axial magnetic field and 
the azimuthal magnetic field, Bq, combine to form a helical screw-type 
magnetic field to confine the plasma. 

Early theoretical studies indicated that the sausage instability 
(m = 0) of a constricted pinch discharge and the kink instability (m = i) 
of short wavelength could be stabilized by an applied external axial 
magnetic field as an alternative to the internal axial magnetic field. 

The possibility of stabilizing a pinched discharge by an external axial 
magnetic field led to the development of che screw dynamic pinch.^ 

The partially stabilized screw pinch can be realized by the appli¬ 
cation of the external axial field simultaneously with the azimuthal 
field of the discharge. This *3 achieved in a linear or toroidal dis¬ 
charge tube of the usual type, modified by imparting a slight helical 
pinch to the return current outside of the tube. The conductor consists 
of a number of copper-braided strips (as in the Tokamak TM-3) insulated 
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from one another, arranged along the length of the tube with a slant 
rdB/dz equal co 0.1-0.2 (i.e,, along the z-axis at radius r and angle 
9). As a result, an external magnetic field is generated by thi 
discharged current, and the total magnetic field (both Bg and con¬ 
tains and heats the plasma in a helical fashion. 

From the facts that the region outside the constricted plasma is a 
partial vacuum and has appreciable conductivity, the concept of screw 
dynamic pinch evolved. To simplify analysis, a pressureleas ionized 
gas of infinite conductivity outside the pinch column proper is assumed. 
The radial equilibrium equation of magnetic pressure for a screw dynamic 
pinch distribution in the plasma column is given by 


dB z b 0 d(xB Q ) 

8 z dr x dx 


(2.7) 


where x = r/r and r = radius of the linear pinch tube.; 
o o 

Given the properties of adiabatic invariant, the magnetic fluxes, 

0^ and 0^, between a given point and the pinch axis in the plasma must 
be conserved. Hence, 0^ has a definite functional dependence on 0^. 
Then, from the orthogonality principle, the azimuthal and axial magnetic 
fields can be expressed in terms of 0g and 0 z (similar to the properties 
of the stream function ana the potential function in fluid mechanics) 

~ 0l ^c/ r o dx ’ B z * ^z^ r o xdx (2.8) 


which satisfv Eq. (2.7) of the pinched plasma. 
Since 


r hy d0 r dx 


(2.9) 
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Introduction of Eqs. (2.8) into Eq. (2.9) yields 

B q = Vt xB (2.10) 

0 o z 

where V - d0^/d0 z which is a constant at any point in the plasma. This 
equation represents the relationship between B^ and B^ at any point x 
in the plasma. In a linear screw pinch, the return conductor of the dis¬ 
charged currents is a helix outside the plasma tube. If the radius of 
the pinch tube is taken as unity, then the wall condition of the screw 
pinch becomes 


( 2 . 11 ) 


r Q x = 1 , 

V Q = (d0 Q /d0 z ) 


* constant. 


x=l 

Substitution of Eq. (2.10) into Eq. (2.7), and rearrangement, give 
the relations 


( 2 . 12 ) 


where 

b - a + v* v>/(i + y 2 0 r V), 

at x = 1, b = 1, and 

B = B 
z zo 

on the wall of the test tube. For given values of r Q and V, the radial 

distribution of B r /B and B /B at any point x of the screw pinch can 
(j zo z zo 


be determined. 
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In the above analysis, the conductivity outside the pinch column has 
been assumed as infinity. Actually, a finite resistivity is present. Con¬ 
sequently, the interdiffusion between Bg and as well as the decay of 
B zq will occur as discussed in the method of internal magnetic field 
stabilization. 

A theoretical study indicated that a screw or diffuse pinch with axi- 

symmetric configuration operating above the Kruskal-Shafranov current limit can 
2 2 

be stable if f) > 0.1 where 3 ■ 2p, Q p/(B z + Bg) and p = plasma kinetic pressure. 
For the purpose of comparison, some computed results of a numerical example 
for the stabilized pinch with external current and the screw pinch with 
external current are shown in Fig. 1. The pinched olasma in both cases 
of the diffuse pinch is stable. The stability is possible for case (a) 
where 3 is in the range 0.05 to 0.10, while for case (b), 3 must be greater 
than 0.1. The magnitudes of the selected magnetic fields B z and Bg are 
typical values and the pitch is defined as P = rlWBg. 

Most of the present laboratory devices with high 3 toroidal plasma 
containment exhibit the flute or interchange instability (m = 11. Dynamic 
stabilization of the m = 1 instability for high 3 screw pinch was studied 
experimentally. 8 In the experimental work, two basic dynamic stabilization 
methods were performed: (a) The oscillating, axially uniform high-frequency 
component of the field, and (b) the oscillating axial component current 
methods. For case (a), special attention is given to the possibility of 
working in resonance with natural plasma oscillation (having relatively 
low frequency) to save energy input required by a high frequency oscillator. 
For case (b) the effects of dynamic shear on the flute instability is 
particularly considered. These methods are similar' to that of dynamic 
stabilization to be discussed later. 





PLASMA TUBE RADIUS x ( - x/x q ) 


Fig. 1 . Typical magnetic field and magnetic pitch of the stabilized 
diffuse pinches vs plasma tube radius: (a) stabilized pinch 
with external current; (b) screw pinch with external current. 


1 
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In the oscillating axial magnetic, field stabilization method, the 
oscillating component of the axial magnetic field sin Ji^t can be pro¬ 
duced and superimposed on the base axial magnetic field B zq in a common 

coil to obtain the total axial magnetic field, B = B + 'B sin a r 

z zo zo st 

where B^ is the amplitude of the oscillating axial field, il is the 
stabilizing frequency and t is the time. The ratio B^/B q can be 
varied between 0.05 and 0.20. The range of uu can be adjusted from 
1-1.7 MHz. The growth rate of the flute instability is moderately 
suppressed after the oscillating axial magne'ic field has been applied. 

In addition, the high frequency energy input could be saved and the plasma 
column breakdown by wall contact could be avoided if the resonance con¬ 
dition between the oscillated plasma motion and the natural radial plasma 
oscillation could, in principle, be achieved. In the oscillating axial 
current, the total axial current, I = I + I sin * t, is fed into 
the screw pinch using the same electrodes for the high frequency oscillating 

I and the quasi-steady current I . T and I have comparable ampli- 
zo J zo zo zo 

tudes from 15-50 kA. The stabilization frequency 0 J st is about 1 MHz. A 

reduction of the m = 1 growth rate is observed after 1 is superimposed 

zo 

on the screw pinch. 

The growth rate or instability frequency, f^ n = uo^ n /2 tt, for the m - 1 
mode instability in a linear (circles) and a toroidal (dots, aspect ratio 
= 5) screw pinch as a function of total axial current \ ? is shown ir 
Fig. 2. The growth rate and the wavelength of the m l instability -are 
very similar in a linear and a toroidal high p screw pinch for adequate, 
corresponding plasma parameters. Neither the end effects nor the elicits 
of toroidal curvature strongly influence the m ® 1 growth rate. The 





GROWTH RATE OR INSTABILITY 



TOTAL ELECTRICAL CURRENT, I , KA 


ig. 2. Growth rate of the m = 1 instability in a linear and toroidal 
screw pinch vs total electric current. 



0.2 0.4 0.6 0.8 l.( 


Typical radial distributions of magnetic fields and electron 
density in a linear screw pinch. 


ELECTRON DENSITY n (10 
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observed radial distributions of the magnetic fields B^, Band the plasma 
density n in a linear screw pinch at t = 2 ^sec, p = 45(10""*) Torr and 
T\ = 0.20 keV is given in Fig. 3. It is seen that the dynamic effect of 
screw pinch can heat the ions of a plasma at very high rates, and the 
dynamic shear tends to stabilize the plasma. The above observed data 

are in good agreement with the calculated results from the theoretical 

. . 9,10 

analysis. 

2.4 Magnetic Shear 

Magnetic shear is appreciably effective in controlling the gross, 
universal, or drift instability of a thermonuclear plasma. It is the 
electrical short-circuit principle whereby any separation of charges from 

positive ions and negative electrons can be neutralized by short-circuiting 

11 12 

the charges. As a result, the plasma can be appreciably stabilized. ’ 

The magnetic shear is produced by the twist of magnetic field lines 

in the plasma. The degree of twist of the magnetic lines can be made to 

vary with distance from the central axis, the farther from the axis, 
the deeper the pitch of the lines, and vice versa. 

The magnetic shear of the magnetic field configuration stabilization 
can be produced by a number of methods. The two most common methods are 
the reversed magnetic field device and the helical magnetic windings. 

The influence of a magnetic shear tends to make the radial pressure 

gradient positive so that Suydam's necessary condition of plasma sta- 
13 14 

bility may be satisfied. * In most case., the magnetic shear is used to 

stabilize the gross, universal or drift instability.^ The gross or 
universal instability (independent of plasma configuration) due to temperature 
and density gradients is synonymous with confinement, and the drift 
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instability is the diamagnetic drift caused by the gradients. 

The drift instability arises from the Landau damping mechanism 
involving particles whose velocities along the magnetic field resonate 
with wave velocity. Sometimes, this damping mechanism also involves 
energy transfer from the resonant particles to the oscillating waves. 

The drift instability can occur under certain conditions:, (a) a large 
density gradient exists in the plasma; (b) a steep transverse particle 
temperature gradient exists; (c) the temperature gradient opposes the 
density gradient of the plasma; or (d) a sufficient, local electric 
field exists or electron current flows along the magnetic field. 

In the introduction of the wave packet model, a local instability 
or perturbation of the plasma is assumed to be composed of wave packets 
$("r, t) of the form 

$(r, t) = 0(r) exp i(E*r - uut) (2.13) 

where 

§Cr) = the spatial perturbation as a function of the space 
vector ?> 

k = the spatial wave number, 

uu = angular frequency of the plasma, and 

t = time. 

A general stability criterion of the wave packets is that the wave energy 
of the wave packets should not grow or exponentiate many times as they move 
through some region of local instability. 

To stabilize the drift instability, the sheared magnetic field in 
relation to the wave packet has two main effects: (a) an effective inrrea 
in the ion Landau damping available for plasma stabilization is through th 
propagation of the' wave packet, or the expansion of the polenti.il well. I 
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potential well can expand into the region where the Landau damping by the 

interaction with resonant ions is strong; and (b) the decrease in wave 

energy of the wave packet, or the shrinking of the wave packet as it 

propagates through some local region and gradually vanishes. 

In stability analysis, a plasma with a density gradient dn/dx 

(n = plasma density) in the x direction, an unperturbed magnetic field 

vector If = B q jlT^ + (x/L^e^j and a Maxwellian velocity distribution is 

considered, where B Q e 2 = B z = main axial magnetic field, B 0 (x/L g )e y 

= shearing magnetic field, Bq = the scalar magnetic field, e y and e z are 

the unit vectors in the y and z directions, and L is the characteristic 
J s 

shear or shearing distance. 

With the Boltzmann distribution function fj(r, v, t) = f^, the Vlasov 
equation is given (in Mks units) by 


df i _ e « - _ _ 

1 + v-Vf. + — 1 (E + v X B)V f = 0 , 
dt j m^ v j 


(2.14' 


where v = velocity vector in space, e = unit electric charge, m = particle 
mass, and E * self-consistent electric field strength. The vector E 
satisfies the Poisson equation of charge (or the Gaussian law) 


• E - L « n 
e o j J J 


in which Gq = permitivity in a vacuum and n^ is evaluated from the jth 


species particle distribution function f 


nj = J^fj(r, v t)d 3 v. 
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If there is no initial equilibrium electric field, E_ = 0, the lineari- 


af 1 + ; • Vf lj + ^ <» * Wlj - - ^ <E i + 7 * V • Voj 


By the method of characteristics, the solution of Eq. (2.15) can be 


J (l i +?X V • Vo. 


where the integral is performed along the particle orbit in the lowest 
order. It is assumed that a time space dependence of the perturbation 


in which (f 1 )q is a constant. By combining the solution of the Vlaso 
equation and the well-known Maxwell equations, the dispersion relation 


for low-frequency modes of the plasma can be given by 
(7 ,-1/2 \ t / ...X t 2 2 


(uu - uuI )Jl 1 - i 


T / uu \ kiiv. 


Is 

ki 1 - yb>e-'’J 4 II B T i 


: electron frequency, 


= electron thermal velocity, 
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: ion frequency. 




ion thermal velocity. 


T e> T. = electron and ion temperatures, or kinetic temperatures, 

2 j \ _ 

k^ = (k^ -! k y ) = wave number perpendicular to B z < = B q 2 z ), 

kjj = k z + j*- k^ = wave number for shear effects, 

b '¥ 2 *l~¥b 2 v 

a. - ion gyroradius or gyromagnetic radius, 
v^ = the Alven velocity, 

l Q (b) = zero order of the first kind of modified Bessel function, 


In deriving the dispersion relation, the following basic assumptions 
have been made: 


• tb: space-time dependence of perturbation varies as the behavior 
of the wave packet defined above, 

• the characteristic length of density gradient, r * = —| 

« k^. The amounts of magnetic shear required to stabilize 
drift instability due only to density gradient for the plasma 
m e^ m i < 0 << ^ is sought. 

• the electron gyromagnetic radius is moch smaller than the wave¬ 
length. 

• ki » k , v. « jiu/k I « v 

i z’ l 1 z 1 e 

• kXp <' 1 and uu « 0^, where X^ = the Debye length and 
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If * T^, 2v^/v^ * 0 (= kinetic pressure/magnetic pressure of the plasma), 
k l >;> k jj » k i ” k2 » and c * k J/ k || = bv A ~ velocit; y of light), then the 
dispersion relation reduces to 



(2.19) 


Based on the dispersion equations, some limiting cases for long and 
short wavelengths of plasma oscillation can be discussed. 



This equation consists of three modes of plasma oscillation: (a) the elec¬ 
trostatic mode, (b) the slow Alfven mode, and (c) the fast Alfven mode. 

2 2 2 

For the electrostatic mode, it is assumed that uj « kyv^ and h 1 • 

The last term in Eq. (2.20) is negligible compared with other terms. An 
obvious solution of the resulting expression is = - i« e as and 
propagate in opposite directions. 

Introduction of iu^ + 6uu^ for a perturbed frequency uu into Eq. (2.20) 

gives 
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Hence, the perturbation condition of tne plasma (as u» e t 0) is 



or 


2ui < kf.v* and 2b ;» kf, vf/u)^ 
o II A ~ I! i e 

which shows that the unstable electrostatic modes propagate almost perpen¬ 
dicularly to the axial magnetic field and appreciable magnetic shear effect 
kjj is required to suppress the drift instability. 

For the continuous (or transient) Alfven modes, the second term on 

, , 2 2 n 

the left-hand side of Eq. (2.20) approaches zero, i.e., u>(u> - u> e " k jj v A " 0< 
The roots of this quadratic equation represent the approximate frequencies 
of the slow and fast Alfven modes: 


Introducing ^ and + 6u> 3 for iu into Eq. (2.20), respectively, the 



(2.23) 


2 2 2 

The perturbation condition for Eqs. (2.22) require ui^ > ^ k j| v A * por a R 1V 
value of v A , a relatively small magnetic shear is needed to inhibit the 
drift instability. The respective frequencies of the electrostatic mode 
and the slow and fast Alfven modes of plasma oscillation are represented 
by u» l + 6a> 1 (= - ou p + 6^), i« 2 F 6 uj 2 and uj 3 + buy 
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2.4b Short Wavelength Oscillation 

In the oscillation of short wavelength b > 1 and lQ(b)e ^ ~ (2rrb) 
Then Eq. (2.19) becomes 


■MS 


l(2nb >" i + fsrfa ■ 1 Ut.,1, 

k fA I!' 


-- 0. (2.24) 


For the Alfven mode, uu^ > 4k|v^ £,n< * t ^ ie ^ as i c assumption « )uu/k|| | <«' v^, 
Eq. (2.24) gives the unstable Alfven mode 

jj / „ uu 2 eAo \!i 


For this case, the effect of magnetic shear kji oecomes insensitive to 
stabilize the drift instability of the plasma. 

2.4c General Shear Stability Criterion for Convective Modes 

To evaluate the spatial growth of a perturbation with frequency ur pro¬ 
duced by a narrow source, the form of the relevant approximate solution is 

expj ik x (x, m)dx 

where k is a wave number and can be found from a dispersion relation, 
x 

This form is equivalent to the expression of the wave packet (Eq. 2.13;. 

When the amplification reached by the perturbation in a plasma region is 

. 19,20 

less than n , the general stability criterion can be expressed as 
x 

max 

- 1m J k x^ X ’ vJu ^ dx < n 0 (2.26) 

o 

in which exp(n Q ) represents a tolerable level of growth of the convective 


(or transient) modes. 
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To illustrate the use of the general stability criterion, Eq. (2.20) 
and the first of Eqs. (2.23) for the convective Alfven mode can be taken 
to advantage. As before, letting uu = uu^ + in Eq. (2.20), it follows 

A 


■ Imk = - 


<]' 


+ 1nTT 

T :»t eg rating this, and considering 6uu 2 as a parameter, the amplification 
of the growth of the propagating wave packet becomes 


k dx = - - 


(\ + * 2 + HV . 
’ \ K "el ' * 


H(f + k W)% + s7ir k ' |V *] 


for simplicity, assigning the value 2u)^ ~ kjj v^, taking the value o 
oob 2 corresponding to the amplification, and performing the integration, the 
stability criterion of the convective Alfven mode is 


where 

L = shearing distance, 
s 

r = - (■— 4^1 = characteristic or scale length of density 

\n dx/ 

gradient (or plasma radius), as defined. 
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which deals with the oscillations and stability of an inhomo&enoons plasma 
contained in helical magnetic field. 




2 2 2 

Recalling the required condition u> e > 4kjjV^, the approximation 
2 

u>„ = iw (kiiv /u> ) and carrying out the integration, the stability 
c e || A e 

criterion becomes 


T-J-.ft)’©''®;)' 


It is important that with high |3 f the magnitude of magnetic shear needed 
for stabilizing the convective Alfven modes is relatively small. 

For purposes of comparison, for higher values of magnetic shear and 
larger values of the relation exists such ^hat XjJuu/kjj vj < 1, in the region 
where the ion Landau damping is strong enc-igh to suppress drift instability. 
In this case, uu - - by a slight modification of Eq. (2.18), the stability 
criterion for the electrostatic modes driven by density gradient is 


which is corresponding to the stability criterion obtained by Krall 
19 

and Rosenbluth 

For the localized region of the perturbation extending into the region 

of the electron Landau damping, Mikhailovskaya and Mikhailovsky obtained, 

21 

however, another stability criterion 


L /r < (r/a V\ 


which deals with the oscillations and stability of an inhomogeneous plasma 


contained in helical magnetic field. 
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1 dn . 
— — and 
n dx 


In the above wave packet analysis, only the density gtadient dn/dx, 
or dn/ndx of the plasma has been considered. If the temperature gradient 

dT/dx, or dT/Tdx (where T = T\ = T g ) is also considered, the parameter, 

dT/Tdx d/dx(An T) d(ln T). . .... .. 

dn/ndJ ‘ d7£(ta7>' or d(ta^> ln the * dircction has “ i " portant cffect 

on the magnetic shear and plasma stability 1 ^. In general, for a large 

ion temperature gradient, no small magnetic shear r/L g is sufficient to 

stabilize the temperature gradient instability where r” 1 

L = characteristic shear distance as defined. 

* 

From the wave packet analysis for both density gradient and tempera¬ 
ture gradient present in the plasma, the following analytical results 
can be obtained. 

• For a density gradient with a nonparallel electron temperature 
gradient d An T g /d An n < 0, or a density gradient with a small 
parallel electron temperature gradient, d An T g /d An n <■ 2, a 
magnetic shear required for plasma stabilization is r/L 

( ' 'i 1/3 

~ (m /m ) 

e i 

• For a density gradient with a large parallel electron tempera¬ 
ture gradient, d An Tj d An n > 2, the magnetic shear required 
for plasma stabilization is r/L ~ (m /m.) 2 . 

SCI 

• For a density gradient with a nonparallel ion temperature gra¬ 
dient d An T^/d An n < 0, or a large parallel temperature gra¬ 
dient d An T^/d An n > 2, there is no magnetic shear related to 
iWm^ ^ oun< * ^ or t * ,e Plasma stabilization. 
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2.5 Magnetic Well 

It has been experimentally demonstrated that universal or interchange 

instability is considerably inhibited in a thermonuclear plasma contained 

22 

in a magnetic well, i.e., at a minimum magnetic field — minimum B. 

Theoretical analysis also predicted that interchange instability of the 
plasma can be contained in a magnetic well, i.e., at a (on the average) 

— 23 24 

well-shaped magnetic field — minimum average magnetic field or minimum B. ’ 
The magnetic well may be conceived as the flux surface of magnetic field 
(or magnetic surface). In other words, the plasma may be contained in a 
magnetic well where the field strength is increased in every direction 
away from the plasma. In general, the magnetic well concerns the shape 
of magnetic field lines. 

2.5a Minimum B (Open-End System) 

The concept of the minimum B is to stabilize a thermonuclear plasma 
in a magnetic mirror device of the open-end system. For the interest of 
a magnetic mirror, if the requirement of self-consistent equilibria and 
plasma potentials in minimum magnetic fields is satisfied, a stability 
criterion for moderate curvature of axially-symmetric magnetic surface 
. 25 

is given by 


- dj l 


1 dB dp 
dx dx 


> 0 


(2.31) 


where p is the kinetic pressure, B is the magnetic field of the plasma, 
and the integration is performed along the element djJ of the magnetic line. 
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In an axially symmetric field, it is adequate to take the magnetic 
flux Y for x as the radial coordinate. Hence Eq. (2.31) becomes 


_1_ dB d£ Q 
4 dY dY 


(2.32) 


Plasma stability is feasible if the pressure is maximum at the minimum B. 
In other words, when dp/dx and dB/dx have opposite signs, Eq. (2.31) or 
(2.32) is positive and the stability criterion is satisfied. 

From plasma energy consideration, a stability criterion related to 
the particle temperatures and the magnetic well depth for univeral 
instability can be obtained. If the initial free energy of the plasma 
is expressed as the Helmholtz function, A, then A(0) is given by 2 ^ 

A(0) - (nVT e )(X 2 /6r) at t = 0 (2.33) 

where 

n = n = electron density, 

V = plasma volume, 

T fi = electron temperature, 

X ■ wavelength, width of local instability zone, and 
r = " (n dx) 33 

bet the instability be localized in a zone of width X equal to a typical 
wavelength. A large number of charged particles diffuse across the field 
lines. If the plasma expands within the zone X(dn/dx), the change in 
average energy is n(dB/dx) (where gi » constant magnetic moment), 

OB = (T, + T e ), and the well depth VB - r(dB/dx), then the initial change 
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in energy, Hv, required for disturbances (in which plasma expands across 
the field in every x direction), is 


As usual, negative 6W implies stability. By combining Eqs. (2.33'' and 
(2.34), the stability criterion for minimum-B to stabilize universal 
instability at the initial state is 


VB _ 1 

B - 6(1 + T./T ) 
i e 


(2.35) 


provided that the oscillation frequency is less than the ion cyclotron 
frequency and the wavelength is greater than the ion gyroradius. If 
T g — T\, the stability criterion for the depth of the magnetic well 
required is reduced to 


VB/B > 1/12 (2.35a) 

if the particle temperatures are not uniform, for simplicity, let 
the particle temperatures T\ ~ xdTj /dx and J’ e ~ nt '"’y point 

x .if the thin plasma column, then the stability criterion !><'< nines 



(2.36) 
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2.5b Minimum Average B (Minimum B, Closed-End System) 

In the toroidal system (dlosed-end system), it is usually found that 
plasma particles diffuse away rapidly across the magnetic field lines as 
in the magnetic mirror system (open-end system). To prevent the universal 
or interchange instability causing plasma particle leakage, the magnetic-well 
idea can be utilized by placing toroidal plasma within a deep magnetic well. 
This idea, however, cannot be done by conduction coils or rods as in the 
magnetic mirror system. 

Although a true toroidal magnetic well is impossible for the closed- 
end system, the net magnetic well or minimum average B (minimum B) for the 
stabilization of a toroidal plasma has been effective.^The trend of 
the minimum average B seems toward more sophisticated and more precisely 
controlled magnetic field configurations. 

A magnetic field uith helical syim/try can provide a mean or average 

magnetic well, i.e., provide regions wherein decreases away from a 

n 2 2 

magnetic axis (or equivalently V = d V/dY is negative where V is the 
volume and 'i' is the flux of a magnetic surface enclosing the axis and the 
volume). The integral can be carried out along the element Al of a magnetic 
field line. In mosc cases, however, the quantity of the second derivative of 
the volume per unit flux, v", is conveniently evaluated in terms of 

v "=If ,= jf <0 < 2 - 37 ' 

z 

in which B^ is the axial magnetic field along the z-axis. 

Stabilization by negative v' (i.e., by a favorable gradient of d//B 
or of average B) is possible in toroidal system and can bo more effective 
than stabilization by magnetic shear against the resistive interchange r.iuieH 
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Minimum average B derived from negative V" of the toroidal system 
against the universal instability is about the same value as minimum 
B of the magnetic mirror system. 

Magnetic wells with central-rod conductor to carry currents (i.e., 
hardcore) may have three types of the stagnation-point (or stationary-point) 
solution: (1) Linear periodic-multipole structure, by using helical coils, 
for instance, 1 = 2 and 4 for shaping field with normal stagnation points, 
and by superimposing 1 = 1 and 3 for corrugating field to create favorable 
VB regions of magnetic well. The basic principle in designing a negative 
V " system is that along any given magnetic field line the favorable curvature 
should outweigh the unfavorable curvature; (2, Helical flux equilibrium, 
by using gross helical curvature to create favorable VB regions and by 
using an axial conductor around which the helical equilibrium flux tube 
is wound; and (3) Toroidal flux equilibrium, by using gross toroidal 
curvature to create favorable VB regions of the magnetic well. For each 

solution, a magnetic scalar potential is properly assumed to satisfy the 

29 

necessary and sufficient conditions of the problem. 

Magnetic wells of the helical field without central-rod conductor 

are also considered. The quantity of the negative V ts evaluated by 

27 

assuming some stream functions for magnetic flux. From a practical 
viewpoint, the central rod is inconvenient. 

2 .b MuLtipolc Configuration 

Multipole configuration consists of two major types linear multi¬ 
pole and toroidal multipole. The linear multiple device has been operated 
as a quadrupole, generally in the collisionless regions, with plasma 
density n = 10*° - 10^ cm ^ and temperature 1-10 eV.- 3 ° Plasma loss is 
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mainly across the magnetic field of the quadrupole. Toroidal multipoles 

have been studied analytically and experimentally from low p equilibrium 

31-TG 

and flute instability to high (3, ballooniig stability. The con- 

9 -3 

finement of low density (n =10 cm ) collisionless plasma with T^ = 40 eV, 
T g = 10 eV produced by gun injection or with T g = 1 eV, T. ^ 1 eV produced 
by microwave heating was performed to determine the mechanisms for plasma 
loss from a toroidal multipole magnetic field. 36 The machine has four 
rods carrying current in one direction and the wall cf the device carrying 
the return current. The plasma loss was attributed to various mechanisms 
of the mechanical hoops, guard supports, local walls, plasma fluctuation, 
diffusion, and the generation of low frequency electric fields by the 
guard obstacles. The chief advantage of the multipole configuration is 
to study plasma behavior in a magnetohydrodynamic (MHD) stable region 
or in the JdA/B stable region. 

To study the plasma confinement, stability, and fluctuations, the 
typical magnetic field shape in quadrupole and octapole configurations 
is shown in Fig. 4. The magnetic field strength distribution in 
the device (along the central vertical line through the median plane and 
along the chamber wall) varying with plasma radius is given in Fig. 5. 

For the uuadrupole configuration, the magnetic line of force which forms 
a cusp, returns around the current-carrying rods. For the octapole con¬ 
figuration, the four rods (with magnetic lines) carrying currents are all 
in the same direction, and the plasma chamber wall carries the return 
current. A stabla plasma region has its maximum pressure on the lines 

t = 0. The last stable or critical line is Y . For stability, the 

c 

plasma pressure must decrease from Y = 0 to i = Y ^. 
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(b) OCTAPOLE 


As pointed out above, two prin¬ 
ciples (or methods) are considered 
important for stabilization: the 
principle of magnetic shear and the 
principle of magnetic well (minimum 
B for the open-end system, and min¬ 
imum B for the closed-end system 
with short connection length between 
the favorable curvature and the 
unfavorable curvature). In the toroi¬ 
dal multiple configurations, the 
magnetic shear can be varied by super¬ 
imposing a toroidal magnetic field. 

The depth of the magnetic well and 
the connection length can be altered 
by changing the number of poles, for 
instance, from octapole to quadrupele. 

The superposition of a toroidal 
field produces both magnetic shear 


* = 0 STABLE PLASMA 

* = CRITICAL STABLE LINE 

-- = MAXIMUM PRESSURE LINE 

Fig. 4. Typical magnetic field 
shape in quadrapole and 
octapole configuration. 


number. The effect of the magnetic 


and finite parallel wave numbers. 

In case the toroidal field B^ is mud' 
smaller than the average multi pole 
field B , the parallel wave number k n 

m II 

is approximately given by kjj - (B^/B^lk^ 
where is the perpendicular wave 
shear tends to stabilize the plasma in 


a multipole configuration. 
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o 



Fig. 5. Magnetic field strength distribution in the device vs plasma 
radius. 

The magnetic well depth of a toroidal multiple configuration can be 
defined as {(Jdjfc/B),^ /2 - (Jdi/B)^. l/( JdA/B),^ where Y is the 
magnetic flux function and Y £ is the flux at the critical line (see Fig. A). 
The magnetic well of the octapole configuration is usually much deeper 
than that of the quadrupole. The flux surface where the magnetic shear 
vanishes is closed to the Y = Y fi surface when the aspect ratio of the 
device is sufficiently large. 

In the toroidal multipole devices, drift cyclotron instability and 
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interchange (or universal) instability have been identified. The 
corporation of the magnetic shear and the magnetic well with short connec¬ 
tion length between favorable and unfavorable regions into the multipole 
configuration may improve the plasma stability. 






3. DYNAMIC STABILIZATION 

The purpose of dynamic stabilization is to perturb the equilibrium 
parameters of a whole thermonuclear plasma to stabilize it. This can be 
accomplished analytically or experimentally by assuming (or applying) an 
oscillating magnetic field, oscillating electric field, or microwave 
heating. The dynamic stabilization is effective, which can stabilize many 
modes of one instability and other types of instabilities simultaneously.^’^ 
The dynamic stabilization of a thermonuclear plasma may consist of 
the following proposed methods: (a) Traveling wave and standing wave, 

(b) oscillating electromagnetic field, (c) oscillating electric field, 

(d) electron ueraperatur" modulation, and (e) axial alternating current. 

Other methods of dynamic stabilization may be developed for certain 
purposes or circumstances once the basic behavior of stabilization 
mechanisms is better understood. 


3.1 Traveling Wave and Standing Wave 


The dynamic stabilization of linear 9 pinch has been previously 
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proposed. In the model under consideration the collisionless plasma 
consists of charged particles which can be reflected into it from its 
surface. Every perturbation of the plasma gives rise to an excessive 
average magnetic pressure which opposes the unstable deformation on the 
plasma surface. In addition to the axial magnetic field, the applica¬ 
tion of a periodic constant azimuthal magnetic field along the pinch 
axis can produce a significant effect on the m = 1 instability of the 
plasma. 
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The dynamic stabilization of the P = 1 bumpy or toroidal 8 pinch 
41 

has also been studied. In the analytical study, the dynamic stabi- 
lizing currents are applied perpendicular to the axial magnetic field. 

(Without the dynamic stabilization, the configuration is basically un¬ 
stable.) The plasma is assumed incompressible because in the analysis 
of the static pinch the most unstable modes of the plasma are found 
incompressible. The new concept introduced in the analysis is to make 
the corrugated magnetic field travel along the pinch axis at the wave 
velocity v^ (where v^ > v^ and v^ = Alfven velocity) because the outward 
drift of the plasma due to the curvature of the torus can be suppressed 
by corrugating the plasma surface. In other words, the corrugated mag¬ 
netic field and surface profile are made to propagate along the pinch 
axis with a wave velocity v w by suitably alternating part of the current 
in the conduction coils. This concept is referred to as "dynamic sta¬ 
bilization of the plasma by traveling wave" (instead of standing or sta¬ 
tionary wave at steady state). It may be assumed that the frequency 
of the magnetic field should be greater than v^/X with the traveling 

wave and less than v /X with the standing wave where v is the sonic velo- 
s s 
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city or the mean ion Velocity and X is the wavelength. Dynamic stabili¬ 
zation by traveling wave is more effective than that by standing wave, in some 
cases. The idea of the dynamic stabilizacion of a thermonuclear plasma 
either by the traveling wave or by the standing wave depends basically 
on the equilibrium principle that all resultant destabilizing forces at 
any point of the plasma (contained in a linear or toroidal system) must 


be vanished. 
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3.2 Oscillating Electromagnetic Field 


The desirability of stabilizing the universal instability of an 

inhomogeneous plasma in a magnetic field by dynamic stabilization has 
40-43 

been studied. The application of an oscillating electromagnetic 

field in stabilizing plasma is briefly analyzed in this section. Con¬ 
sider a collisionless plasma confined in the x direction by a uniform 
magnetic field B z> parallel along the z axis. The scale length of 
density gradient is given by r * = - dn/ndx where n is the plasma density 
(n decreasing with increasing x). The temperature gradient of the plasma 
is small conpared to the density gradient. An oscillating electromag¬ 
netic field which consists of the electric field E and the magnetic 
field B can be expressed in terms of the scalar potential 0 and the vector 
potential A 

E = - 70 - * ^ B = V X A (3.1) 


where 


k 


y’ 


0 O (X), 


A = a^e* a 1 (x) + ...), 

0 = a o (e i6 0 1 (x) + ...) + a g0 Q (x) , 

0 = k^y + k z z + lit, 

0 = the frequency of the electromagnetic field, 
c = the speed of light, 

k = wave numbers along the y and z axes, respectively, 

(x) = given functions, and 

a = amplitude of the electromagnetic field and is assumed 
to be small, and a^x) are given vectorial functions. 
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The stabilizing effect produced by the electromagnetic field is due to 
the drift velocity of the guiding center of ions and electrons. The 
dynamic stabilization is effective when the drift velocity is of the 
order of the diamagnetic ion velocity v rf . (= negative diamagnefc elec¬ 
tron velocity v de , where - T./m.™., T. = T = ion temperature, 
m. = ion mass and uu ci » ion cyclotron frequency). 

In the absence of the oscillating electromagnetic field let o 

’ c 

be the radius of gyration, v the velocity parallel to the z-axis, 
and x £ the abscissa of the guiding center of a charged particle. The 
motion of the particle is not in resonance with the electromagnetic 
field. 

In the presence of the electromagnetic field, the motion of the 
charged particle at a given point x is specified as 


V + W V a " d V + Vv V 


T = (Q + k. V + k v )t, 

2 y y z. z' ’ 

V *z 3nd *c " the mean velocities in the y and z directions and the mean 
cyclotron frequency of the particle, respectively, 
t = time, and 

f( T l» T > = the periodic function. 

For given time t and particle energy W, the mean velocities \ the mean 

cyclotron frequency of the particle as functions of p , v n x can be 

c II c 

given by 
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V y<P C > V || ’ 

V-.- v ll- 

“c <p c’ v ll • 


* > 


1 



dU 

dx 


c 


x ) 


1 dw 

m bv|| 
_2_ dW 


M) * 2 

c a P c 


(3.2) 


qB z 

where i« c = —, q = charge, and m = mass of the particle. 

For an inhomogeneous plasma with uniform temperature T, the scalar 

potential 0 T may be written as 

i(k y+k z4u>t) 

0j. = b(x, 0)e y Z 


in which b(x, 0) is a function of x and 0, periodic nature in 0 with a 
period 2TT. Comparing the expressions given for Eq. (3.1), b(x, 0) can 
be expanded into series as 

b(x, 0) = b Q (x: + b 1 (x)e i9 + b^e'’ 19 + ... - b Q (x) 


where b^x), b ^(x), ... are small compared with b Q (x) and are of the 
order of magnitude corresponding to the small amplitude, a . 

ik x 

For drift modes of low frequency « |^ ci |» b Q (x) = h(x)e x 

where h(x) is defined by the condition Jdh/hdx| < l/4r^ (where 
2 \ 

r. = [2T./m.uT.jp = mean radius of gyration of ions), the dispersion 
i l i ci 

relation of the plasma oscillations in the presence of the oscillating 
44 

electromagnetic field is given by 

L. + L + Dal = 0 (3.3) 

l e 0 

where L(L^ and L g for ions and electrons) is related to the mean value 


of the particle 






36 


+ Vd>’o<Vc> M‘^u >h 2 (x c )d ^n dx c 


f Q (p c ,v^ ) or f(p c »V|j ) - the Maxwellian distribution exp m(uu^p^ + )/2 t| 


J (k.p ) = the Bessel function of zero order of the first 
ole 

kind, and 

D = real and bounded in absolute value independent 
of the modes. 

The plasma is dynamically stabilized if the real parts ReL^ > 0 and 

ReL e > 0 when |Da^| « 1, |uu| < |kyVj| , and |dh/hdx| < l/4r^. These 

conditions can be satisfied if the values of function v (p , v , , x ) 

y c il c 

for each set in the domain, p 2 < 2T/mu) 2 , v 2 < T/m and lx I < r. are 

H c c’ || 1 c> 1 

spread over an interval of the order of || ( = = diangnetic velocity). 

c 

There are two cases related to the dynamic stabilization by the 
standing wave of the oscillating electromagnetic field: 

(a) when the oscillating electromagnetic field is an electromagnetic 
standing wave along the x-axis, then 

k = k = 0, lu) I « Q» k (T/m ) z and U» w (uj = (4nnq 2 /'m ) ? 
y z 1 c 1 xe pp ^e 

= the plasma frequency). 


electromagnetic field induced by an exterral source, it follows that 


0 (x) 
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Hence Eqs. (3.2) are reduced to 


v (p , v , , x ) 
y V || c 


k 212 

s _x_ q e 

muu „ n 2 ‘'o'^’W 
c 2m 12 


J (2k p ) sin 2k x 


where E is the mean square intensity of the electric field of the oscillating 
electromegnetic field. 

(b) when the oscillating electromagnetic field is an electrostatic 
standing wave along the x-axis: as before, let k z = a^(x) 

= 0 o (x) = 0 

|(2/u) c jJ « |k^/k x J « 1, and 0^(x) = 2b cos k^x 

If the electromagnetic field is produced in the plasma, then 
Eqs. (3.2) become 


W v v - ■ 


—^ J^(k p ) cos 2k x 


3.3 Oscillating Electric Field 


Oscillating high-frequency electric field applied for dynamic sta¬ 
bilization of an inhomogeneous plasma consists of two cases: 

(a) Uniform high-frequency electric field, when the skin laver 

6^ of the field is equal or greater than the transverse dimensions of 
the plasma L^, i.e., when 6^ > L^, the electric field is assumed 
uniform. 

(b) Non-uniform high-frequency electric field, when the skin layer 
depth < L^, the electric field is considered non-uniform in 
the plasma. 

In both cases the stabilizing effects depend on the plasma parameters 
and the frequency and strength of the electric field. 
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A thermonuclear plasma confined by a magnetic field in a closed-end 
device is subject to both hydromagnetic and drift instabilities. These 
instabilities can result in anomalous diffusion of the plasma across the 
magnetic field. High-frequency dynamic stabilization of hydromagnetic 
instability has been studied^ dynamic stabilization of drift (gross 

or universal) instability by oscillating electric field has also been 

. . . 46-48 

considered. 

3.3a Uniform High-Frequency Electric Field 

The stabilization of collisionless drift instability in an inhomo¬ 
geneous plasma is feasible by the application of a uniform high-frequency 
electric field E(t) in the direction parallel to the magnetic field, 
which may be called the "E-stabilization": 

E(t) = E cos Ot (3.4) 

where E(0) = E q , t is the time and Q is the frequency of the applied 
electric field. When a high-frequency electric field of sufficient 
amplitude operates in the plasma, the oscillations of electrons in the 
field become important. If a drift instability is excited in the plasma, 
the plasma electrons are subject to an additional force in the direction 
parallel to the magnetic field due to the presence of the high-frequency 
electric field, E(t) and E^(tj, *)• This force is equal to the electron 
charge times the mean values of the electric field, e(E + E^), where E^ 
is the part of high-frequency induced electric field. The parameter 
dUn T)/d(£n n)plays a very important role in the determination of sta¬ 
bilizing effect for different frequency regions related to the parameter 
of ion speed ratio, (here T = plasma temperature and n = plasma density 


as defined). 
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When this additional force is in phase with the force due to kinetic 
pressure of the plasma, the electric field in the drift wave is increased 

and so is the frequency, uu. Under the conditions (3\/nu) 2 « uu/k^ « (T^/n^) 2 
k 2 2 

and = (k^/uu L )(2T i /m i ) of the intermediate-frequency case (kj^ = k^ 

+ ky, = Larmor frequency), as uu is increased, the Landau damping 
by the electrons is intensified. Therefore the collisionless drift 
instability of the plasma tends to be stabilized. The range of values 
of the parameter d(jlnT)/dUn n)in which high-frequency stabilization 
occurs is maximal for -• 0 and approaches zero when -♦ ». For the 
low-freqeuncy case, uu « k (T^/nij) the instability region is the 
same as in the absence of high-frequency electric field d (in T)/d(Xn n) 

> 2(1 + 6) where 6 = ^[l - I x <C?/2)/I Q (C?/2)] and l o <pJ/2), l^pJ/2) 
are the modified Bessel functions. 


j <•' k Q C i (T i /m i )^ and k z « k^^/m^ are given and the requirement of t lie- 
parameter d(^nT)/d(jtnn) < (1/C ± - 1) is satisfied, the dynamic stabilization 
of a collisionless inhomogeneous plasma by an oscillating uniform high-frequency 
electric field with modest amplitude can be achieved. The value of k Q is 
defined as k Q = (n o T)" l d(n Q T)/dx, where n Q and T are the plasma density and 
kinetic temperature at equilibrium state. 

The electron-ion collisions in a collisional, inhomogeneous plasma 

resulting in the excitation of drift instability is known as the drift- 
49 

dissipative instability. When a dispersion relation is derived on the 

basis of the usual hydrodynamic equations, the collision term in ( , v 0 ^ 

due to electron-ion collisions must be taken into account in the equation 

of electron motion parallel to B (where m is the electron mass, v is 
1 z e e 

the - electron velocity and v is the electron-ion collision frequency in 
the 7. direction).. It can be shown from solutions of the dispersion 
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equation that the drift-dissipative instability can also be stabilized 
by a high-frequency electric field. The amplitude of the uniform high- 

frequency electric field required, however, is higher than that for 

47 

the collisionless drift instability obtained. 

3.3b Non-Uniform High-Frequency Electric Field 

As discussed, if the skin layer depth 6^ of the high-frequency 
electric field is comparable with the transverse dimension L A or plasma 

radius, the high frequency electric field can be considered as uniform 

47 

and the effect of the magnetic field can be neglected . If o A «- L A 
the high-frequency field skin effect, however, must be taken into account 
The high frequency electric field is then non-uniform and the effect 
of the magnetic field component becomes important. In such case, the 
stabilizing effect of the non-uniform high-frequency electric field 

48 

gives rise to two additional independent mechanisms of stabilization : 

(a) the gravitational mechanism produced by the averaged high-frequency 
field forces in an unperturbed plasma state, (b) the oscillating 
magnetic-field mechanism of the time-dependent component of the wave 
field vector projected in the axial magnetic field B^ direction. 

The gravitational mechanism of stabilization may be called the "f-stabili- 
zation," and the oscillating magnetic field mechanism of stabilization 
nay be called the "B-stabilization." Thus, the effects of the B-stabili- 
zation, E-stabilization and g-stabilization of drift waves by a non-uniform 
high-frequency electric field can occur concurrently on a weekly inhomo¬ 
geneous plasma confined in a constant, main magnetic field, B^. 
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Consider a collisionless, inhomogeneous plasma confined in a strong 
constant magnetic field B^ q along the z-axis. The applied non-uniform 
high-frequency electric field can resolve into two components: 


= E (x) sin fit, and 
o 

„ dE (x) 


(3.5) 


where E q (x) is the amplitude of the applied electric field E, is the 

induced magnetic field component parallel to the main magnetic field 

B , and fi, t and c are as defined before. The plasma density n and 
zo 

temperature T vary with x radially. 

The averaged high-frequency field force which influences the 
g-stabilization in an unperturbed plasma state is given by 

F «V-HEr <3 - 6) 


where 


g = the gravitational acceleration, 

Y = the high-frequency potential function, 

E = E (x) in the z direction, and 
zo o 

a. = e, i for electrons and ions in the plasma respectively. 
The high-frequency field forces act mainly with the negative electrons. 
In the steady state, however, there is charge separation together with 
a potential electric field E = - V0 (see Eqs. 3.1), or = - o0/dx. 

As a result, the positive ions are also contained. 

The conditions of a collisionless, inhomogeneous plasma require 
that the dimensions of particle density gradient, r 1 = - “ and 
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particle temperature gradient, (r 1 > T = - ~ and the skin layer depth 6^ of 
the high-frequency electric field, are large compared to the Larmor radius 
of the ions. 

From the oscillating electric and magnetic fields and the field force 
produced by the non-uniform high-frequency electric field and the con¬ 
ditions required for a collisionless, inhomogeneous plasma confined in 
a constant magnetic field, the dispersion relation can be derived and 
solutions of the dispersion equation can be obtained. In the absence of a 
high-frequency electric field, the plasma is stable only against a fast 
ion acoustic wave when the parameter d(£nT)/d(ton) > 0. The plasma is 
unstable for d(£nT)/d(£nn) < 0. When the high-frequency electric field 
is applied, the stable region broadens toward negative values of 
d(£nT)/d(£n n). The effect of stabilization is the increased frequency 
Ui which accompanies the applied high-frequency electric field, as a 
result, the Landau electron damping is intensified. Thus, the stabilizing 

effect produced by the high-frequency field pressure i? associated with the 

8 a 

Doppler shift due to drift velocity v, = -* (where ui = qB /m ) of the 

Q U>^ Cx Z 0( 

charged particles, i.e., the mechanism of the g-stabilization. 

In general, the stability of a collisionless, inhomogeneous plasma 
stabilized by a non-uniform high-frequency electric field depends on 
the plasma parameters and the frequency and strength of the electric field. 

An effective stabilization of the plasma by a non-uniform electric field 
may be achieved near the resonance condition. 

Analytical work to predict the effects of dynamic stabilization on 
drift instability of an inhomogeneous plasma by an oscillating electric 
field parallel to the magnetic field has been discussed above. Very 





43 


recently, experimental work to observe the effect of dynamic stabili¬ 
zation of the plasma by the oscillating electric field has also been 
reported.In the experimental work reported, dischage-type helium 
plasma and Q-machine cesium plasma have been used. Plasma diagnostics 
are made by means of radially movable electrostatic (or Langmuir), and 
magnetic probes. Electron density, electron temperature, magnetic field 
strength, modulation frequency, wave amplitude, ac voltage, etc., are 
measured. The decreases of drift wave amplitude and wave spectrum in 
several modes (m = l, m = 2, ...) prove the stabilizing effect of the 
electric field at relatively low and high frequency regions. 

A proposed schematic diagram for dynamic stabilization of a dis¬ 
charged plasma by a high-frequency electric field is shown in Fig. 6. 
Basically, the system consists of two cathodes (with a tubular anode 
between them) in a magnetic field parallel to the axis of the anode 


QUARTZ 

A NOW TUBE 



ELECTRIC FIELD 

Fig. 6 . Schematic diagram for dynamic stabilization of a dischargedplasm 
by high-frequency electric field. 
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tube. Deuterium gas is discharged through perforated tantalum cathode 
at the right, and plasma emerges from a tabular cathode at the left. 

Ions produced in the discharge are accelerated back and forth between 
the two cathodes, where the ions produce secondary electrons. These 
electrons are prevented by the magnetic field from reaching the anode, 
and they in turn produce more ions in the plasma. Ne-'tral particles are removed 
by the divertor. A sinusoidal ac potential with variable amplitude and 
frequency is applied near the left end of the plasma column. The average 
current of the plasma is measured by ammeter A, Fig. 6. Plasma diagnostics 
can be made by means of movable electrostatic probe, magnetic probe, 
spectroscopic measurement, visual analyzer and microwave techniques. 

3.4 Electron Temperature Modulation 

The effect of a small, periodic, time-dependent modulation of the 
electron temperature on the drift waves of a collisionless, inhomogeneous 
plasma is considered. The energy source required for electron temperature 
modulation can be a periodic, intensity-varying microwave to heat the 
plasma. If the transit time of a particle through the plasma device is 
short compared to the modulation period of microwave heating, the elec¬ 
tron temperature will follow the modulation of the energy source. Energy 
transfer from the electrons to the ions in the plasma will be small if 
the electron-ion equipartition time is long compared to the modulation 
period. Ohterwise, the energy transfer from the electrons to the ions 
will be large. 
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If the plasma has a Maxwellian velocity distribution, the modulation 
electron temperature, T^x, t), may be expressed as 

T (x, t) = T (1 + a cos fi t) (3.7) 

e eo m 

where T is the peak electron temperature, a is the parameter a « 1, 
eo 

and is the modulation frequency of the energy source such as micro- 
wave heating. By using the Vlasov equation of electron motion and 

giving a small perturbation of the equilibrium parameter, a dispersion 
52 

relation can be derived. For the plasma with weak density gradient, 
temperature gradient and/or axial current in the limit of .anishing Larmor 

radius, the stability criterion for different modes of drift wave can be 

. 52,53 
given by 





46 


Equation (3.8) is valid as long as a « 1 and v i « |j^-| « v . It is 

z 

seen that the plasma parameters, d(in T eo )/d(jfcn n)and dOfcnTjJ/dUn n) 
remain to play the important role to stabilize the drift waves by electron 
temperature modulation. 

The theoretical results show that the dynamic stabilization of 
drift waves by electron temperature modulation for a collisionless 

inhomogeneous plasma is possible if the stability criterion E^. (3.8) 

* 

is satisfied and one of the beat frequency, is of the 

order of |kjv^ where the ion Landau damping becomes intensified. 

3.5 Axial Alternating Current 

Dynamic stabilization of collision-dominated drift waves by means 
of an axial alternating current (or ac current) is studied. With the 
collision-dominated region, an inhomogeneous plasma has a uniform mag¬ 
netic field B(= B^) in the z direction and the density n and the temperature 
T varying in the x direction. When an axial electric voltage is applied 
to the grid of the plasma tube in the presence of an ac current, the 
drift waves of the plasma at low frequencies can be suppressed. 

Drift waves in a cesium plasma are stabilized experimentally by 

applying an external, ac current parallel to magnetic fields through 
54 

a meshed grid immersed in the plasma. The results of the dynamic 
stabilization shown in the experimental data may be interpreted by the 
facts: (a) The bias effect due to rectifying characteristics of a 

plasma sheath when ac voltage or current is applied to the grid, 
causes the electron current to increase because of rectifying charac¬ 
teristics of the electron-rich plasma sheath, and (b) some non-linear 




47 


Landau damping effect due to couplings between drift waves and external 
perturbations of the applied current when the applied current frequency 
0 is closed to drift wave frequency, and strong coupling between them 
is observed. 

The simple, analytical result for dynamic stabilization of collision- 
dominated drift waves in the presence of an ac current is also obtained 
when the applied current frequency Q is comparable with drift wave 
frequency.^ A macroscopic electron velocity u q (x, t) parallel to the 
magnetic field can be produced by the ac current flowing in the 
plasma. From the equation of electron motion, the linearized equation 
of continuity, the energy balance equation of the system, and 
simplification, the stability condition for low-frequency drift waves is 
given by 


>U o d T eo Uv 2 

, 2d£nn o kV 

eo L z e< 


.2 2 d £n n 

k v o 

z eo 


/ d Xn T \ / 


1 T - \ 
_io)> 

1 n o/ 


(3.9) 


2*i a i 


v = (2 T /m )\ 
e e e ’ 

V = the electron-ion collision frequency, v ’ and 

o = the initial or unperturbed state of the plasma.; 

It is seen from Eq. (3.9) that the combined effects of finite ion 
Larmor radius, the temperature gradient and density gradient of ions and 
electrons, the currer- gradient, and the electron inertial influence are 


all involved in the stability condition of collision-dominated drift waves.. 
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4. FEEDBACK STABILIZATION 

Plasma stabilization by feedback methods has been introduced quite 
recently.'’*’ **** The selected and proposed feedback methods to stabilize a 
thermonuclear plasma are: (a) electrostatic feedback; (b) neutral atom 
beam feedback; (c) modulated microwave energy sources, and (d) modulated 
uniformly distributed electron sources. Some new concepts of oscillating 
or helical magnetic fields for feedback stabilization of a thermonuclear 
plasma can be developed similarly. 

4.1 Electrostatic Feedback 

Feedback stabilization consists of sensing the plasma instability at 
any point in the plasma by sensor devices and then feeding back an 
accurate signal of proper amplitude and phase in order to suppress the 
instability. A number of feedback stabilization experiments have been 
conducted with electrostatic (or Langmuir) probes as sensors as well as 
control elements. Such experiments have proved successful in suppressing 
collisional drift instability of the plasma. Ihe theory of electrostati" 

feedback stabilization of drift waves to interpret the experimental 

, . , , . J 56-60 

results has also been advanced. 

Various electrostatic feedback t chniques have been developed for 
plasma stabilization. By modulating a parallel-electron-current sink, 
for example, (i.e., modulating and drawing electron current from an 
electrostatic probe parallel to the axial magnetic field) the collisional 
drift wave is stabilized with negative feedback in the unstable region 
of the potassium plasma at T = 0.23 ev, n Q = (10^)cm \ dn^/n^dx = - 1.25 cm 
and the critical magnetic field B c = 2.19 kilogauss (kG). The drift 
instability is the m = 2 azimuthal mode. The experimental results obtained 
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for constant feedback gain varying with phase (angle) shift of the 
applied probe (suppressor) voltage for various magnetic fields are shown 
in Fig. 7(a-c). Figure 7(a) shows the plasma density or instability 
amplitude, Fig. 7(b) gives the change of plasma density at the center 
of the plasma column, and Fig. 7(c' chows the instability frequency as 
functions of feedback phase shift in both stable and unstable regions 
respectively. The theoretical instability growth rate and frequency 
of the plasma are also shown in the curves of Fig. 7(a) and (c). The 
feedback gain varying with the plasma density, characteristic length of 
the density gradient and instability frequency of the plasma is shown 
in Fig. 8(a-c). The radial distribution of equilibrium plasma density 
n Q , plasma density or instability amplitude n, and required minimum 
probe (suppressor) voltage V for feedback stabilization is given in 
Fig. 9. The plasma density at the center of the plasma column as a 
function of magnetic field strength with and without feedback stabilization 
of collisionai drift waves is shown in Fig. 10. 

With a suitable feedback circuit arrangement, stabilization of the 

m = 2 aximuthal mode can be achieved up to B = 1.2B (B = critical magnetic 

c c 

field strength). Thus, the electrostatic feedback (negative or positive) 
demonstrates that collisionai drift instability driven by density gradient 
can be controlled from penetration of feedback signals into tne plasma 
interior with phase shift. 

4.2 Neutral Atom Beam Feedback 

The application of feedback stabilization by means of electrostatic 
P’-obe or material suppressor is limited to plasmas of moderate density 
,.nd temperature because an electrode cannot resist very high temperature 
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PLASMA RADIUS, mm 

Radial distribution of plasma density, n 0 , instability amplitude 
(density), n, and required minimum probe (suppressor, rms), vol¬ 
tage, V, for feedback stabilization. 
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Plasma density at the center of the plasma column .as a functior 
of magnetic field B with and without feedback stabilization of 
collisional drift waves. 
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of the plasma. For a thermonuclear plasma, the injection of a neutral atom 
beam to provide feedback-controlled volume sources of particles and momen¬ 
tum density is a basically different approach from the electrostatic feed¬ 
back. Any perturbation imparted to the plasma by a neutral atom beam can 
be sensed by optical or microwave techniques. 

By using the equations of motion of ions and electrons (i.e., two-fluid 
equations for the plasma), and the equation of continuity with source terms 
of density and momentum, after linearization and simplification, a quadratic 
equation of the dispersion relation for frequency id is given by^ 

k.v v 

uu(uu - UJ. - llh.) + itD (id - U) - i(D- + -*—nD-) + —= 0 (4.1) 

it s 6 t uj t r r 

c 

where 



T. 

v di = = ion diamagnetic drift velocity, 

T 

V de = qBr = electron diamagnetic drift velocity, 

uu^ = S/n^ = feedback frequency or gain, 

v = neutral atom beam velocity, 
x 


M = neutral atom mass, 
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= equilibrium plasma density, 

= neutral atom density, 

S = neutral atom beam source, 

t) = plasma resistivity related to oj 

s 

T\,T e ,q,B,k ,k z> etc., have been defined. 

For o; s » |uu e | and = 0, Eq. (4.1) reduces to 
2 

uu - du(uu. - iuj ) - iuu UJ =0. 
l s s e 

This has the two solutions 

w = uu + iu; (uj - u).)/u) , and uu, = jj. - uu - iu> (4.2) 

ae e e is bie s 

The purpose of the analysis is to use the neutral atom beam to stabilize 

without destabilizing ua^. 

There are three basically different stabilization mechanisms which 

are characterized by the phase (angle) shift 9 and the relationships between 

v and ru (or normalized beam momentum v /riti ) of the neutral atom beam 
x e x e 

feedback system:' 

• Plasma density smoothing, i.e., smoothing density perturbation 
with the neutral-beam injected plasma. For this case 

)v^| « |tOJ e | and 0 = 180°. The marginal stability requirement 

in Eqs. (4.1) and (4.2) is US- = ui (uu - x) )/uu. For x, » lx I 
r s e s 1 e 1 

» lw I, this becomes iu, = m (uj - uu.)/ju = - v = the instability 
1 f 1 f e e i s T 

growth rate. 

• Equivalent (or simulating) minimum B configuration |v^| » ( r j^| 
and h = 0°. The marginal stability condition is = - J c ^' c » " * j * 
The pressure exerted by the momentum of the neutral beam for this 





case becomes important and produces an equivalent or simulating 
minimum B effect. 

• De-energization of resistive drift waves. For this case 

|v^| » |raJ e |j 0 = - 90°. The marginal stability requirement 
is the feedback gain = iyo^/kjv^. The momentum of the 
neutral atom beam exerts pressure which reduces the energy of 
the waves in the plasma. 

Comparison of the effectiveness of the three mechanisms indicates that 
feedback gain required for stabilization by mechanism (a), density smoothing, 
is rather large. Momentum control of the neutral beam, mechanism (b), can be 
more effective than density control if v^ — v^. At the same time, mechanism 
(b), as equivalent minimum B configuration may have a relatively long 
energy replacement time of the plasma. The alternative momentum control, 
mechanism (c), can be better than mechanism (b) for very small instability 
growth rate, y/u^, hut it is less effective than mechanism (a). 

1 2 

In a numerical example, it is assumed that a> e = - > 0, b = —(k^a^) 

= k 2 . —- = 2(10 -2 ) or 10" 6 , 0 = arg(ai c /uj ), uu /uu = 20 or 200, v /r* = 650 
•12 °fese xe 

MJj 

c 

- normalize beam momentum. By introducing these values into Eqs. (4.1) 
and (4.2), the calculated values of instability growth rate <md minimum 
feedback gain required for stabilization are obtained. Figure 11 shows 
the growth rate varying with phase shift for various values of the nor¬ 
malized feedback gain 0^/111^. Figure 12 indicates the minimum feedback gain 
I"f/^e^min re 9 uirec * f° r stabilization varying with the normalized neutral 

beam momentum v /run . It is seen that the value of |e r /j I decreases 
x e 1 1 e 1 m 1 n 

as . / 1 . increases when the feedback stabilization for resistive drift 
s e 


wave is concerned. 
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4.3 Modulated Microwave Energy Sources 

The possibility of feedback stabilization for resistive drift waves 
of a thermonuclear plasma by the remote injection of neutral atom beams 
has been analyzed above. In this section some experiments to demonstrate 
remote feedback control of collisional drift waves by modulated microwave 
energy sources is discussed. The experiment is performed in a cesium plasma 
of the Q-machine with T g = 0.24 eV, = 5(10^) cm ^ and B = 4.0 kG. 

A half-wave dipole antenna (located about 2 cm outside the plasn' olumn) 
is used as a remote suppressor, and a Langmuir probe is use as .. detector. 
The amplified and phase-shifted ion-saturation current instability signal 
is introduced to modulate the microwave energy source so that the micro- 
wave power output is proportional to the instability amplitude. The 
microwave energy heats the plasma by resonance absorption. 

For optimum stabilization, the modulated microwave power must be 
applied to the region of maximum instability amplitude and approximately 
180° out of phase with the local density oscillation. The measured 
results from the plasma for (a) the dependence of instability amplitude 
and frequency on feedback phase shift and gain, (b) the radial distri¬ 
bution of plasma density, and (c) the change of Magnetic field strength 
on plasma stabilization with microwave feedback are similar to those 
obtained by the electrostatic feedback (see Figs. 7-10). 
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4.4 Modulated Uniformly-Distributed Electron Sources 


An inhomogeneous plasma confined by a main magnetic field B z in the 
axial direction is considered. The particle density n varies in the x- 
direction. The feedback stabilization with modulated electron sources uniformly 
distributed in the plasma is applied to a very weakly collisional regime 
where the electron mean free path is greater than the parallel wavelength. 

In this case the motion of electrons can be described by fhe kinetic equa¬ 
tion in the drift approximation. The kinetic equation for perturbed electron 
distributions with collision and source terms in the simplest form is 
given by 


E df 
h _y _oe 

B dx 


df 

;_01 

dv 


df. 
_le 

dz 


= - Vj£.~ n 


where 


f = f( v", t) = the Boltzmann distribution function (see 
section 2.4) 

f , f. = the zero-order and first-order distribution 
oe le - 

functions 


E = electric field in the y-direction 

y 

m, q = the electron mass . charge 

, n. , v , v = the electron densities, velocity, and collision 
e 10 e frequency 


f = f /n = the Maxwelliam velocity distribution. 
M oe oe 


The feedback electron source, S , is distributed throughout the plnsra and 
e 

responds linearly to the local electron density perturbation. Based on 

, . . . , 59,60,69 

this assumption, the source term can be given by 


(4.4) 
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By the Fourier component analysis of Eq. (4.3) with the aid of Eq. (4.4), 
solving for f. g and integrating, the ratio of the electron density at local 
perturbed state to the electron density at initial equilibrium state. 


n. /n , is obtained 


i - t|, e ) J - Vz + M 

T 1 - i(» £ - vj j {V^/to - + ii^jj 


) = q0 = work function, 0 = scalar potential, 


T = kinetic temperature, 
k T 


k , k , u) = have been defined previously. 


In the case of short mean free path of electrons (for relatively strong 


n - <"* + i(k v _ ) /v. 

le _ _ z Te ei 


where v^ g = (2 T/m) , electron thermal velocity, and v gi = electron-ion 
collision frequency. In the case of long mean free path of electrons 
(for relatively weak collisions), v « k v and oj « k v , Eq. (4.5) 


1 ¥ i(ir/2)^ (a) - uj* + iv ) /(Ik lv ) 


If the ion density is also perturbed by the electron sources, the rati 
of ion density at local perturbed state to ion density at initial eqni 


li _ uu* - txt) 
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1 2 2 

where b = - k a. as defined previously. 

For b z: “>£/ k z v Te « 1, from the quasineutrality condition the dis¬ 
persion relation is obtained. 



(4.9) 


Feedback stabilization can be achieved with purely imaginary feed¬ 
back frequency uu^ such that arg(u) f /uo*) = - 90° (or + 270°), and gain 
required for 

|uj f u)*| > 2b (id*) 2 (4.10) 

Feedback stabilization can also be achieved with phase shift 0° if the 
condition 



is satisfied. Thereofore, the optimum phase shift for the feedback sta¬ 
bilization will be between - 90° (or + 270°) and 0°. The electrostatic 
feedback by phase shift of - 90° and the neutral atom beam feedback by 
phase shift of 0 are corresponding to the phase shift of the modulated 
electron source feedback. Although the approaches of the feedback stabili¬ 
zation are different, in general, the feedback stabilization by the mod¬ 
ulated electron sources is in good agreement theoretically and experi¬ 
mentally with the corresponding stabilizing effects by the electrostatic, 
neutral atom beam momentum or modulated microwave energy source feedback. 
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5. MAIN COMPARISONS AND CONCLUSIONS 


5.1 Main Comparisons 

In the magnetic field configuration stabilization, the internal 
axial magnetic field, the rotational transform, and the external axial 
magnetic field of screw (bumpy or diffuse) pinch together with the 
azimuthal magnetic field to form some fundamental elements of magnetic 
fields for plasma stabilization. The magnetic shear and the magnetic 
well (including minimum B for open-end system and minimum average B for 
closed-end system) have shown good results to control the hydromagnetic 
(MHD) and drift instabilities. The multipole configuration consists of 
the linear and toroidal multipoles. The linear multipole device has been 
operated as a quadrupole in the collisioniess regime. The toroidal mul¬ 
tipoles have been studied from low (3 equilibrium and flute instability 
to high p-ballooning stability. The magnetic shear, the magnetic well, 
and the multipoles can be considered as the advanced magnetic field 
configuration for plasma stabilization. 

The high magnetic power, p mag » required per unit length of the confinement 
chamber for these configuration is approximately given by 


P 

mag 


~ C 


b 2 rti 

NA 


(5.1) 


B = the magnetic induction, or magnetic field strength, 

” = the coil radius, 

A = the cross-sectional area of the coil, 

N = the number of coil turns per unit length of the confinement 


chamber. 






7] = the resistivity of the coils, and 
C = a constant. 

Among the dynamic stabilization methods all use oscillating elec¬ 
tric field or alternating current, except the electron temperature modu¬ 
lation method, to stabilize the drift waves of a thermonuclear plasma. 

In the oscillating electromagnetic field method and the oscillating non- 
uniform electric field method, both electric and magnetic fields are 
considered. In the electron temperature modulation method, 
microwave (or radiofrequency field) and optical energy sources can be 
employed to heat the plasma by remote control. The oscillating magnetic 
field for use in the dynamic stabilization is usually produced by alter¬ 
nating current or currents. 

Of the methods discussed above, ‘•he dynamic stabilization for drift 
instability of a collisionless, inhomogeneous plasma is most effective. 

The axial alternating current method used for dynamic stabilization, how¬ 
ever, is potentially effective for collision-dominated drift waves of 
the plasma. 

Dynamic stabilization has demonstrated the possibility of stabilizing 
low-1requency drift waves (or broad turbulences) and improving contain¬ 
ment time of a thermonuclear plasma. Apart from electron temperature 
modulation by microwave heating, modulated neutral beam or optical energy 
source (such as weak laser beam) may be useful for dynanic stabilization 
of a fusion reactor by the remote control and detection methods. 

The power required for dynamic stabilization, by the high- 
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frequency electric field at frequency U is given by 
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vhere 

n Q = equilibrium plasma density, 

T, b = as defined previously, and 

bn Q T = a fraction of the interrial (or thermal) energy of 
the plasma. 

Electrostatic feedback for which feedback elements immersed in, or 
in contact with, plasma have proved effective to stabilize plasma insta¬ 
bilities and improve plasma confinement. Particularly, the collisional 
drift instability in the oscillatory regime has been stabilized by immersing 
electrostatic probes to draw electric current from the plasma. Since 
material probes cannot resist high temperature in a thermonuclear plasma, 
the electrostatic feedback method is probably limited to laboratory 
plasmas of moderate density and temperature. 

The application of feedback control to a nuclear fusion reactor, how¬ 
ever, depends critically on the possibility of remote control and detection 
method. Therefore, the remote control of neutral atom beam feedback to 
stabilize resistive drift waves and the modulated microwave energy sources 
to stabilize collisional drift waves has been proposed. Charged particle- 
beams and optical beams (including laser beam) can also in- used for remote 
control and detection in the feedback stabilization processes. 

The effectiveness of some experimental results on the low-frequency 
dritt waves and the good qualitative agreement with theoretical expecta¬ 
tion indicate that the study of feedback stabilization (as well as dynamic 
stabilization) may be a fruitful approach to increase the containment 
time of a thermonuclear plasma. 
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The power required for linear feedback stabilization, P f ^> by the 
parallel electron current is given by 

P„ = a 2 bn Tuu. (5.3) 

fb o f 

in which uu^ = uu^ = instability frequency, a = q0/T ~ n/n = normalized 
instability amplitude, and bn Q T remains a fraction of the internal (or 
thermal) energy of the plasma. 

The ratio of power required for the feedback stabilization to that 
required for the dynamic stabilization on the basis of Eqs. (5.2) and 
(5.3) is 

w ■ *V n - (5 - 4 > 

which is much smaller than unity when a « 1 and < ft in the most 
practical cases. 

The high magnetic power requirement of the advanced magnetic Lie Id 
stabilization system is equivalent to the total magnet loss per unit 
length of the confinement (or vacuum) chamber due to heating, confinement, 
and stabilization of the plasma. 

For purposes of comparison and conclusions, the main advantages and 
disadvantages of the plasma stabilization systems are given below: 

5.1a Magnetic Field Configuration Stabilization 

• Provision for plasma heating, confinement and stabilization 

• Self-stabilization in appropriate magnetic field configuration 
» Suppression of the hydromagnetic (or rtllD) instabilities in 


particular 
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• Complex stabilizing magnetic fields required 

• Complex winding and high cost of the solenoidal conductor 

• High power required for magnetic field generation and operation 

5.1h Dynamic Stabilization 

• Possibility of stabilizing many modes, different low-frequency 
drift instabilities, and broad-based turbulences simultaneously 

• Suppression of low-frequency drift instabilites in particular 

• Stabilizing frequency much higher than instability frequency 
and no phase shift needed 

• Excessive power required to maintain new dynamically stable 
equilibrium for high growth-rate instabilities. 

5.1c Feedback Stabilization 

• Relatively simple magnetic geometries required when appropriate 
feedback and phase shift systems are selected. 

• Lower power requirement for stabilization because when the 
amplitude of instability is reduced, the feedback amplitude 
reduces accordingly. 

• Feedback can interact only with one unstable mode near marginal 
instability; it is difficult to stabilize many modes simultaneously. 

• Application of feedback control by means of material probes and 
electrodes limited only to plasmas of moderate density and 
temperature. Remote control and detection methods are preferred. 

• Linear feedback (feedback frequency U)^ = instability frequency 
with specific phase angle) can stabilize different drift insta¬ 
bilities and the m = 1 flute instability W 
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Nonlinear feedback (repetition rate < instability frequency, with 
specific phase angle) effective to suppress perturbations in 
the nonlinear feedback constant corrective signal is applied 
whenever the instability amplitude grows above a prescribed 
level. Higher modes, however, are easily destabilized. 


Conclusions 

In the magnetic field configurations, the principles of magnetic 
shear and magnetic well will continue to play an important role 
in plasma stabilization. The trend of the magnetic shear and 
the magnetic well is toward more sophisticated and more precisely 
controlled economic magnetic field geometry to stabilize hydro- 
magnetic and drift instabilities. 

Dynamic stabilization (fi » u> in , no phase requirement) can 
stabilize many modes, different low-frequency drift instabilities, 
and broad-band turbulences simultaneously. Excessive power, how¬ 
ever, is required to maintain new dynamically stable equilibrium 
for high growth-rate instabilities. 

Feedback stabilization (selected frequency and specific phase shift 
required) can interact only with a part of the plasma and one 
unstable mode near marginal instability. Feedback control by 
means of material probes and electrodes is limited to plasmas 
of moderate density and temperature. Linear feedback (u> f = »» in ) 
can reduce drift instabilities effectively and suppress the 
m = 1 mode of flute instability reasonably. Nonlinear 
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feedback (repetition rate < instability frequency is 

effective to suppress low-frequency drift instabilty at repetition 
rates below instability frequency but higher modes are easily 
destabilized. 

• There are more basic differences between dynamic stabilization and 
feedback stabilization. In the dynamic stabilization the high 
frequency (or radio-frequency) field interacts with the entire 
plasma at a new equilibrium density n Q while in the feedback 
stabilization the plasma equilibrium remains unchanged and the 
time-varying suppression signal interacts only with the unstable 
part of the plasma. Both dynamic stabilisation and feedback sta¬ 
bilization, however, will prefer to have remote control and 
detection methods for use in a controlled fusion reactor, even a 
rather advanced and costly technique is required for remote 
dynamic or feedback control. 

• At present, plasma self-stabilization in appropriate economic 
magnetic field configurations is still uncertain. If the plasma 
self-stabilization can prove to be fully effective in the labora¬ 
tory scale, the actual design of a controlled fusion reactor with 
the provision of complex stabilizing magnetic fields will re¬ 
quire some reliable scaling laws to determine the reactor size. 

The reliable scaling laws for controlled fusion reactor design, 
however, have not satisfactorily been found. Since the magnetic 
field configuration stabilization is mostly effective in sta¬ 
bilizing the hydromagnetic instability. The dynamic or feed¬ 
back stabilization which can stabilize the drift (.gross or 
universal) instability of a thermonuclear plasma may contribute 
to controlled fusion research. 
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• It is expected that the research and development of remote con¬ 
trol and detection method for potential dynamic stabilization 
or feedback stabilization of a thermonuclear plasma will pro¬ 
bably become important in developing a controlled fusion reactor. 
Either the dynamic stabilization or the feedback stabilization 
may be combined with the magnetic field configuration to sta¬ 
bilize a thermonuclear plasma so that the minimum requirement 
for a controlled fusion reactor may be attained and satisfied. 

In the preceding discussion and analysis a number of new ideas which 
can be developed into original methods for understanding plasma instability 
and improving plasma stabilization have aeen introduced. The main fields 
of the new ideas may consist of the combination of magnetic shear and 
magnetic well, and the combination of magnetic field configuration and 
dynamic or feedback stabilization with remote control and detection methods. 
In the main fields, some specific ideas such as the combined effect of 
rotational transform,magnetic shear and minimum B on a toroidal plasma, 
extention of the present magnetic shear analysis, dynamic stabilization 
by remote microwave (cr radiofrequency) field control (see Fig. 6), etc., 
can be developed to make contributions for a controlled fusion reactor. 
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